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We study the phase behavior of colloidal suspensions the solvents of which are considered to be binary liquid 
mixtures undergoing phase segregation. We focus on the thermodynamic region close to the critical point 
of the accompanying miscibility gap. There, due to the colloidal particles acting as cavities in the critical 
medium, the spatial confinements of the critical fluctuations of the corresponding order parameter result 
in the effective, so-called critical Casimir forces between the colloids. Employing an approach in terms of 
effective, one-component colloidal systems, we explore the possibility of phase coexistence between two phases 
of colloidal suspensions, one being rich and the other being poor in colloidal particles. The reliability of this 
effective approach is discussed. 
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I. INTRODUCTION 

In colloidal suspensions the dissolved particles, typi- 
cally micrometer-sized, interact via effective interactions 
which arc a combination of direct interactions, such as 
electrostatic or van der Waals interactions, and indirect, 
effective ones due to the presence of other, smaller so- 
lute particles or mediated by the solvent^ For example, 
by adding non-adsorbing polymers to the solution deple- 
tion interactions between colloidal particles can be in- 
duced which are predominantly attractive^ The range 
and the strength of this entropy-driven attraction can 
be manipulated by varying the polymer-colloid size ratio 
or the polymer concentration. The resulting phase dia- 
grams of the colloids are very sensitive to changes in the 
depletion-induced pair potential between colloids It is 
also possible to induce such attractive depletion forces 
by adding much smaller colloidal particles or by using 
solvents with surfactants which form micelles acting as 
depletion agents^ 

In contrast to depletion forces, solvent mediated in- 
teractions between colloidal particles can depend sensi- 
tively on the thermodynamic state of the solvent £ This is 
the case if the solvent exhibits fluctuations on large spa- 
tial scales such as the fluctuations near the surfaces of 
the colloidal particles associated with wetting phenom- 
ena near a first-order phase transition of the solvent or 
the thermal fluctuations of the solvent order parameter 
near a second-order phase transition. In the latter case, 
critical fluctuations are correlated over distances propor- 
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tional to the correlation length £ which diverges upon 
approaching the critical point of the unconfined solvent 
(without colloids) at the critical temperature Tc & \ Ex- 
posing the near-critical fluid to boundaries, e.g., by in- 
serting colloidal particles acting as cavities, perturbs the 
fluctuating order parameter near the surfaces of the col- 
loids on the scale of £ and restricts the spectrum of its 
thermal fluctuations. Since these restrictions depend on 
the spatial configuration of the colloids, they resuh— in 
an effective force between the particles which is called 
the critical Casimir force (CCF) fc- Accordingly the 
range of the CCF is proportional to the bulk correla- 
tion length £. Therefore it can be tuned continuously by 
small changes of the temperature T. The range can also 
be controlled by varying the conjugate ordering field 
of the order parameter of the solvent such as the chemi- 
cal potential in the case of a simple fluid or the chemical 
potential difference of the two species forming a binary 
liquid mixture. The strength and the sign of fc can be 
manipulated as well. This can be achieved by varying the 
temperature T or hb and by suitable surface treatments, 
respectively Compared with other effective forces be- 
tween colloid particles CCFs offer two distinct features. 
First, due to the universality of critical phenomena, to 
a large extent CCFs do not depend on the microscopic 
details of the system. Second, whereas adding depletion 
agents or ions changes effective forces irreversibly, the 
tuning of fc via T is fully and easily reversible. 

Although the CCFs between two colloidal particles de- 
pend on the (instantaneous) spatial configuration of all 
colloids, one can consider dilute suspensions or tempera- 

(s) 

tures sufficiently far away from T c , such that the range 
2i? + £ of the critical Casimir interaction between the col- 
loids of radius R is much smaller than the mean distance 
between them. For these cases the assumption of pairwisc 



2 



additive CCFs is expected to be valid. Using this pair- 
wise approximation we are able to map the actual system 
of a mixture of (monodispersc) colloidal particles and sol- 
vent molecules to an effective, one- component system of 
colloidal particles, in which the presence and influence 
of the other components of the mixture enter through 
the parameters of the effective pair potential. Adopt- 
ing this approach allows us to use standard liquid state 
theory in order to determine the structure of an ensem- 
ble of colloidal particles immersed in a near-critical sol- 
vent and to study its sensitivity to changes of the critical 
Casimir potential due to temperature variations. Since 
within this effective approach the feedback of the colloids 
on the solvent and its critical behavior is neglected, this 
way the phase behavior of the full many-component sys- 
tem cannot be determined in all details. However, for 
those parameters of the thermodynamic phase space for 
which this approach is applicable, we predict a colloidal 
"liquid" -"gas" phase coexistence, i.e., the coexistence of 
two phases which differ with respect to their colloidal 
number densities. 

The necessary input for the approach employed in 
the present study is the CCF fc between two spheri- 
cal particles. It is known that at the bulk critical point 
fc is long-ranged. In the so-called protein limit corre- 
sponding to D/R, £/D ;§> 1, where R is the colloid ra- 
dius and D the surface-to-surface distance between the 
colloids, the so-called small-sphere expansion^ renders 
f c (D;T^ s) ,R)/ (fc£>T c (s) ) - R d-2+r, D -(d-i+r,)^ where d 

is the spatial dimension and r\ is the standard bulk crit- 
ical exponent for the two-point correlation function. In 
the opposite, so-called Derjaguin limit D <C R one has^ 
fc{D-M s) ,R)l UbT^A ~ R (d-i)/2 D -(d+i)/2_ Thus 

the CCF can indeed successfully compete with direct 
dispersion^ or electrostatic forces in determining the sta- 
bility and phase behavior of colloidal systems. Away 
from the critical point, according to finite-size scaling 
theory (see, e.g., Ref. [lj), the CCF exhibits scaling de- 
scribed by a universal scaling function which is deter- 
mined solely by the so-called universality class (UC) of 
the phase transition occurring in the bulk, the geome- 
try, and the surface universality classes of the confin- 
ing surfaces.— ~— The relevant UC for the present study 
is the Ising UC with symmetry-breaking boundary con- 
ditions. For spherical particles, theoretical predictions 
for the universal scaling function $ of the CCF in the 
full range of parameters are available only within mean- 
field theory (MFT)ii^ (for h b = and ellipsoidal par- 
ticles see Ref. [Til ). In contrast, for planar surfaces in 
d = 3 results beyond MFT are available. For a van- 
ishing bulk ordering field hj, = and symmetry break- 
ing surface fields they are provided by field-theoretical 
studies^ the extended de Gennes-Fisher local-functional 
method,—^ and Monte Carlo (MC) simulations ,22i22r2& 
Moreover, in Ref . results for the scaling function of the 
CCF for hb ^ are presented which are based on a den- 
sity functional approach. Also corresponding experimen- 



tal data£~— i 28 i 29 are available. Based on the Derjaguin 
approximation 3 ^ the knowledge of the scaling function of 
the CCF acting between two parallel surfaces can be used 
to obtain the scaling function for fc between two spheres 
and between a sphere and a planar walh 17 ' 18 The Der- 
jaguin approximation assumes that the surface-to-surface 
distance D between two spheres is much smaller than 
their radius R. In many cases, however, this approxi- 
mation works surprisingly well 8 -^ even for D < R. Wc 
shall use this approximation for temperatures which cor- 
respond to £ < R, because under this condition the CCFs 
between the colloids act only at surface-to-surface dis- 
tances D small compared with R. In order to handle 
the dependence of the CCF on hb within the Derjaguin 
approximation, wc propose a suitable approximation for 
the film scaling function of the CCF in d — 3. A neces- 
sary input for this latter approximation is the mean-field 
scaling function for the CCF, which we have calculated 
using the field-theoretical approach. 

Experimental studies of the phase behavior of colloidal 
suspensions with phase separating solvents have encom- 
passed silica spheres immersed in water-2-butoxyethanol 
(CaE\ ) mixtures and in water-lutidine mixturesi 32 ' 33 In 
Ref. l3J both silica and polystyrene particles immersed in 
these binary mixtures have been studied focusing on the 
formation of colloidal crystals and its relation to aggrega- 
tion phenomena. A few theoretical 3 ^ and simulatio n 36 ' 37 
attempts have been concerned with such kinds of col- 
loidal suspensions. 

Our paper is organized such that in Sect.[II]we summa- 
rize the theoretical background of our analysis. In Sub- 
sect. Ill Al we discuss colloidal suspensions and the effec- 
tive model we use, whereas Subsect. Ill Bl provides infor- 
mation concerning critical phenomena and the CCF. For 
the parameters entering into the effective potential be- 
tween the colloids, in Subsect. IIII Al we discuss the range 
of their values corresponding to possible experimental re- 
alizations. In Subsect. IIII Bl the thermodynamics of the 
considered colloidal suspensions is analyzed. A general 
discussion of the phase diagram of the actual ternary 
mixtures consisting of the colloids and the binary solvent 
is given in Subsects. IIII B II and IIII B 21 The reliability 
and the limitations of the effective, one-component ap- 
proach are discussed in Subsect. IIII B 31 The results for 
the phase diagrams emerging from the effective approach 
are described in Subsect. [ill B 41 In Sect. [IV] we conclude 
with a summary. 

II. THEORY 

A. Colloidal suspensions 

1. Interactions 

The effective interactions between colloidal particles 
are rich, subtle, and specific due to the diversity of mate- 
rials and solvents which can be used. Our goal is to pro- 
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vide a general view of the effects a critical solvent has on 
dissolved colloids due to the emerging universal CCFs. 
Therefore we adopt a background interaction potential 

(s) 

between the colloids which is present also away from T c 
and captures only the essential features of a stable sus- 
pension on the relevant, i.e., mesoscopic, length scale. 
These features are the hard core repulsion for center-to- 
centcr distances r < 2R and a soft, repulsive contribution 



V re p (r) /(k B T) = U rep (r) 



Aexp(-nD) , D 



2R> 0, 



(1) 



which prevents coagulation favored by effectively attrac- 
tive dispersion forces. The main mechanisms providing 
U rep (r) are either electrostatic or steric repulsion, which 
are both described by the generic functional form given 
by Eq. ([T)) i 38 i 39 The steric repulsion is achieved by a poly- 
mer coverage of the colloidal surface. If two such covered 
colloids come close to each other the polymer layers over- 
lap, which leads to a decrease in their configurational en- 
tropy and thus to an effective repulsion. Concerning the 
electrostatic repulsion, for large values of the surface-to- 
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surface distance D the effective interaction between the 
corresponding electrical double-layers at the colloid sur- 
faces dominates; this leads to a repulsion. The range k _1 
of the repulsion is associated with the Debye screening 
length in the case of electrostatic repulsion and with the 
polymer length in the case of the steric repulsion. The 
strength A of the repulsion depends on the colloidal sur- 
face charge density and on the polymer density, respec- 
tively. For the effective Coulomb interaction screened by 
counterions A is given by 3 ^ 



A={ee a y 1 T 



R/ i.k B T) , 



(2) 



where e is the permittivity of the solvent relative to vac- 
uum, eo is the permittivity of the vacuum, T is the sur- 
face charge density of the colloid, and k is the inverse 
Debye screening length. 

In our present study we shall analyze the behavior of 
monodisperse colloids immersed in a near-critical solvent. 
We treat the solvent in an effective way, i.e., we do not 
consider the full many component system but an effec- 
tive one-component system of colloids for which the pres- 
ence of the solvent enters via the effective pair potential 
V (r = D + 2R) /(k B T) = U (r = D + 2R) (see Eq. Q): 



U(r) 



Urep H~ Ui 



(d=3) 



D < 

Aexp(-KD) + (1/A)6>( d = 3 )(6, A, £), D > 0. 



(3) 



In Eq. ([3]) C/i^ -3 ' is the critical Casimir potential and 
$( d =3) j g j-(- s un i versa l scaling function in spatial dimen- 
sion d = 3. For the present system, 0W is a function of 
the three scaling variables O = sgn(t) D/£, A = D/R, 
and £ = sgn{h b )£/^ h \ Here f(*^0) = fjf 1*1"*, 

with t = ±(T - T c (s) )/ r i S) for an upper (+) and a 
lower (— ) critical point, respectively, is the true corre- 
lation length governing the exponential decay of the sol- 
vent bulk OP correlation function for t — > ± and for 
the ordering field, conjugated to the OP, hb = 0. The 
amplitudes £^ (± referring to the sign of t) arc non- 
universal but their ratio £o~/£cT ^ s universal. The corre- 
lation length ^C 1 ) = \hb\^ v ^^ S ^ governs the exponen- 
tial decay of the solvent bulk OP correlation function for 
t = and hb —> 0, where l;^ is a non-universal ampli- 
tude. (Note that £q has not the dimension of a length, 
see Appendix [A] v, (3, and 6 are standard bulk critical 
exponents^) The critical finite-size scaling and the scal- 
ing variables will be discussed in detail in Subscct . IIIBI 

We do not consider an additional interaction which 
would account for effectively attractive dispersion forces. 
Effectively, dispersion forces can be switched off by using 
index-matched colloidal suspensions. As will be discussed 
in Sect. IIII1 the presence of attractive dispersion forces 
does not change the conclusions of the present study. For 



I 

a particular experimental realization A, K, and R are 
material dependent constants. 



As mentioned in the Introduction, the scaling function 
9^ of the critical Casimir potential between two spheres 
is not known beyond MFT for the full range of the scaling 
variables O, A, and E. In order to overcome this restric- 
tion we shall use two approximations for 9^ d=3 \ First, 
we use the Derjaguin approximation (c.f., Eq. (|23p) in 
order to express 9^ in terms of the universal scaling 
function of the CCF for the film geometry. The lat- 
ter is known rather accurately from MC simulations in 
d = 3 and for = 0; we shall use this knowledge in 
the present study. Second, in order to be able to capture 
the dependence of CCFs on hi, we propose an approxi- 
mation for i?y as a function of S. Our approximation 
is constructed in such a way that for hb — > the scal- 
ing function reduces exactly to i?f| (S = 0) for all d, and 

at Tc^ its shape is the same as within MFT (sec, c.f., 
Eq. ([M]) ). These two approximations will be discussed 
extensively in Subsect. Ill Bl 
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2. Thermodynamics and stability 

Determining the thermodynamic properties of a sys- 
tem from the underlying pair potential V (r) of its con- 
stituents is a central issue of statistical physics.— In prin- 
ciple the thermodynamic properties of a system can be 
determined from its correlation functions . 41 ^ 42 For exam- 
ple, the so-called virial equation provides the pressure p 
of the homogeneous system: 

2 f°° 

p/(pk B T) = 1 - -TTp / U' (r) g (r) r 3 dr. (4) 
Jo 

The radial distribution function g (r) is related to the to- 
tal correlation function (TCF) h (r) according to h (r) = 
g (r) — 1. The isothermal compressibility xt follows from 
the sum rule 41 - 

lim S(q) = pk B Txr, (5) 

q— S-0 

where the structure factor 

S(q) = l + ph(q) = l/(l-pc(q)) (6) 

can be determined by scattering experiments. The 
system becomes unstable for xt — ► oo, correspond- 
ing to the critical point and, within MFT, to the 
spinodals in the phase diagram (see, e.g., Ref. 0). 
The TCF h (r) is related to the direct correlation 
function (DCF) c (r) and the number density p via 
the Orstcin-Zernicke equatio m 41 i 44 h (r) = c (r) + 
pjh (r') c(|r — r'|) d 3 r'. The correlation functions can 
be calculated iterativcly . 41 ' 45 For a given approximate 
expression for a (r) one obtains hi (r) in Fourier space: 

h(q) = c(q)/(l-pc(q)), (7) 

with h (q) ~ Je iqr /i(r) d d r, analogously for c(q), and 
q = |q|. By choosing a bridge function b(r), the closure 

h(r) + 1 =exp{-U{r) + h(r) -c(r)+b{r)} (8) 

renders a DCF c^+i which typically differs from c;. This 
procedure is continued until satisfactory convergence is 
achieved. The initial guess Ci=i is guided by the shape of 
the direct interaction potential. For the bridge function 
we use the so-called Percus-Yevick approximation (PY), 

b PY (r) = In [h (r) - c (r) + 1] - h (r) + c (r) , (9) 

and the so-called hypernetted-chain approximation 
(HNC), 

b HN c(r)=0. (10) 

In terms of the DCF c (r) and of k (r) — h (r) — c (r) 
(which is more useful than h (r) for handling the hard 
core U(r < 2R) = oo) the PY closure (Eq. ([9])) can be 
written as 

cpy(r) = [exp{-[/(r)}-l](fc(r)+l) (11) 



whereas the HNC closure (Eq. (JTOj)) leads to 

chnc (r) = exp {-U (r) + k (r)} - fc (r) - 1. (12) 

For a more detailed discussion of the applicability and re- 
liability of this integral equation approach (IEA) we refer 
to Ref. |42|. For comparison with MC simulations in the 
case of a potential with attractive and repulsive parts see 
for example Ref. H(| which discusses particles interacting 
with a pair potential containing attractive and repulsive 
Yukawa-like contributions Eie~ KiV jr. The IEA is capable 
to reveal the rich phase behavior of such systems. 

There are further relations expressing thermodynamic 
quantities in terms of correlation functions which are ex- 
act from the formal point of view. However, because the 
bridge function is not known exactly, the set of equations 
([7]) and ([5]) may have no solution in the full one-phase re- 
gion of the thermodynamic phase space; moreover the re- 
sulting thermodynamic quantities depend on the scheme 
taken. This is the well known thermodynamic inconsis- 
tency of this approach (although there are more sophisti- 
cated schemes trying to cope with this problem) 42- More- 
over, for the kind of systems considered here, the deter- 
mination of phase equilibria is even more subtle because 
due to the adsorption phenomena, which are state depen- 
dent, the effective potential between the colloids depends 
on the thermodynamic state itself. Inter alia this implies 
that the effective potential acting between the particles 
should be different in coexisting phases. This feature 
is not captured by the effective potential approach pre- 
sented above. Therefore within the IEA only certain es- 
timates for the coexistence curve can be obtained. For a 
reliable phase diagram actually the full many component 
mixture (such as the binary solvent plus the colloidal 
particles) has to be considered. 

Another useful first insight into the collective behav- 
ior of attractive (spherical) particles is provided by the 
second virial coefficient 41 - 

poo 

B 2 = 2n (1 -exp{-U {r)})r 2 dr. (13) 
Jo 

Beyond the ideal gas contribution it determines the lead- 
ing non-trivial term in the expansion of the pressure 
P (p) I (ksTp) = 1 + B2P + • • • in terms of powers of 
the number density p. Measurements of B2 for colloids 
immersed in near-critical solvents have been reported in 
Ref. H3. Vliegenthart and Lekkerkerker— and Noro and 
Frenkel 4 ^ (VLNF) proposed an extended law of corre- 
sponding states according to which the value of the re- 
duced second virial coefficient = ^/-S^ 5 ' 1 at the 
critical point is the same for all systems composed of 
particles with short-ranged attractions, regardless of the 
details of these interactions. B 1 ^ 8 ^ = ^f-cr 3 is the sec- 
ond virial coefficient of a suitable reference system of 
hard spheres (HS) with diameter a . This (approximate) 
empirical rule is supported by experimental data 4 ^ and 
by theoretical results^ The critical value B^ c can be 
obtained in particular from the Baxter model for adhe- 
sive hard spheres^ for which the interaction is given by 
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cxp {— U (r)} = H (r — a) + -^5 (r — a) where H (r) is 
the Heaviside function and S(r) is the delta function. The 
reduced second virial coefficient is related to the so-called 
stickiness parameter S by 



/;. 



2jt 



1 



(14) 



This model exhibits a liquid-vapor phase transition as 
function of S with the critical valued S c ~ 0.113, so that 
B* c ~ -1.212. 



As discussed in Subscct. lIII Al undcr typical experimen- 
tal conditions the bulk correlation length £ and therefore 
the range of the CCFs are smaller than the radius R of 
the (micron-sized) colloids. Therefore the main attrac- 
tive contribution to the resulting effective pair potential 
between the colloidal particles is localized at a range of 
distances D which arc small compared with R. Therefore 
for experimentally realizable conditions the effective pair 
potential can be considered to be short-ranged and thus 
matches the conditions described above for the approx- 
imate VLNF conjecture. Accordingly, close to the cor- 
responding value of £>2 c the effectively one-component 
system of colloidal particles can be expected to exhibit a 
critical point terminating a "liquid" -"gas" phase separa- 
tion line. 

In order to obtain a suitable HS reference system, fol- 
lowing Weeks, Chandler, and Anderse n 53 ' 54 we split the 
pair potential U, as it is commonly done,— into a purely 
attractive contribution 



Ua(r) 



U (r) , r min < r, 



(15) 



where U has its minimum at r„ 



and into an effective 



HS core, exp{ — Ujjs Ml — H (r — a), with the diameter 
<7 defined by 



= / (1 - cxp{-J7(r)}) dr, 
'o 



(16) 



where U (r = r$) = 0. We have checked that for the 
potential U considered here (c.f., Eq. ([26]) ) and for its 
split potential 



U s (r) = U HS (r) + U a (r) , 



(17) 



as given by Eqs. (|T5|) and ([16 
are almost the same. 



the resulting values of Bi 



3. Density functional theory 



Density functional theory (DFT) is based on the fact 
that there is a one-to-one correspondence between the 
local equilibrium number density p(r) of a fluid and a 
spatially varying external potential acting on it. It fol- 
lows that there exists a unique functional J-[p] which is 
minimized by the equilibrium one-particle number den- 
sity and that the equilibrium free energy of the sys- 
tem is equal to the minimal value of the functionali 5 - 5 - 
Since for static properties the absolute size of the par- 
ticles does not matter, DFT has turned out to be very 
successful not only in describing simple fluids but also 
colloidal suspensions^ While the ideal gas contribution 
to the functional J- is known exactly, J- id [p] /(fc^T) = 
J dr p (r) {in [\ 3 p (r)] — l}, where A is the thermal 
wavelength, the expression for the excess contribution 
J~ ex = J~ — J~id is known only approximately. 

As a first step, for a liquid in a volume v we consider 
the simple functional 



J^ex [p] = 
Fhs ([p] ;ct) 



d 3 nd 3 r 2 pCriJC^triaMra), 

(18) 



where r\2 = |ri — 1*2 1 . For the attractive contribution 
U a to the interaction potential entering into Eq. (|18p wc 
employ Eq. (fl5|) . Ths {[p]',&) is the excess functional 
for the HS system for which we consider the effective HS 
diameter a as given by Eq. (Tf6|) . For Fhs various sophis- 
ticated functionals are availablei 5 -^— For our purposes, 
however, it is sufficient to determine the bulk free en- 
ergy density. According to Eq. (fT8]) the so-called random 
phase approximation for J- is given by 



7TCr 3 /6 



Frpa/ (kBT) = T\ a In 



1 - Vo 



2 - IOt?^ + 5t£ \ 1 
2(1-^) 2 ) 2 



(19) 



where for the HS-contribution in Eq. (|T9|) the PY- 
approximation (as obtained via the compressibility equa- 
tion) has been adopted^ U aQ — -^U a (q = 0) with 
U a (q) as the Fourier transform of the potential. The 



r 



packing fraction r] a used in Eq. f| 19|) corresponds to the ef- 
fective HS system, i.e., r\„ = ^cr 3 p = (jji) V- I n Eq. ([I 



Fid/ (k B T) = vp In (n„) - In § (<r/A) J - 1 is incor 



6 



porated, except for the term up In (cr/A) 3 j which ac- 
counts only for a shift of J- linear in p and therefore 
is irrelevant for determining phase coexistence. For the 
free energy given in Eq. (|19p the critical point is given 

implicitly by r, ( a R c PA) = 0.129 and (u a ^ = -21.3. By 

choosing the PY-approximation we are able to compare 
the results obtained by DFT and by the IEA on the same 
footing. 



studied in detail for the film geometry in d — 3 by MC 
simulations) 25 ' 26 in d = 2 by using exact solutions ; 60 ' 61 
and in d = 4 within Landau-Ginzburg theory^ In a sys- 
tematic perturbation theory in terms of e = 4 — d (with 
d* — 4 as the upper critical dimension for the bulk Ising 
UC) thermal fluctuations are captured with a statistical 
wcighiii 3 - ~ exp (— H [(/>]). MFT corresponds to the low- 
est order (d = 4) in this expansion and accordingly the 
MFT-cquilibrium configuration 4>mft minimizes R [(f)]. 



B. Critical phenomena in confined geometries 



The critical Casimir force 



In this subsection, within the field-theoretical frame- 
work we provide the relevant theoretical background for 
critical phenomena in confined geometries. We pay spe- 
cial attention to the case of a binary liquid mixture close 
to its demixing point. 



1. General concept 

The free energy J 7 of a system (i.e., the solvent in the 
case considered here, see below) close to its critical point 

\Tc S \c^j is the sum of a regular, analytic background 



contribution J> and a singular part J- s . Within the ap- 
proach of the field-theoretical renormalization group the- 
ory the leading behavior of the singular contribution !F S 
for a confined system is captured by the (dimensionless) 
effective Landau-Ginzburg Hamiltonian 



n[(t>] 



hb4> 



d d r 



(20) 



where v is the volume available to the critical medium 
and dv is its confining (d — l)-dimcnsional surface. <f> — 
4> (r) is the local order parameter (OP) field associated 
with the phase transition. For the case of binary liquid 
mixtures, which in their bulk belong to the Ising UC, the 

OP 6 /c a s ,c is a scalar. The quartic term 

with the coupling constant g > stabilizes the Hamil- 
tonian for > r oc t (t = ± {t/T [ c s) - l) , sec the text 

below Eq. ([3])) and hb is a symmetry breaking 6ulk field 
conjugate to <j>. In Eq. (f20f the surface couplings en- 
compass the so-called surface enhancement c(r) and the 
surface field h s (r) with r G dv. Equation (f20|) turns out 
to capture the fixed-point Hamiltonian of surface critical 
phenomenal For laterally inhomogeneous substrates h s 
and c vary along dv£L For colloids strongly preferring 
one of the two species of the binary liquid mixture the so- 
called strong adsorption limit applies, which is described 
by the so-called normal fixed point, i.e., (h s = oo, c = 0) 
for all r € dv. Corrections to the fixed-point behavior of 
the CCF due to finite surface fields and the crossover 
between various surface universality classes have been 



For the case of colloidal suspensions considered here, 
the colloids act as cavities in the critical medium and 
thus in Eq. ([20]) the confining surface dv is the union of 
the surfaces of all colloids in the system. Accordingly, 
the OP profile for a given colloid configuration depends 
on the position of all colloids and therefore the CCFs 
which act on the confining surface dv are non-additive. 
In order to cope with this very demanding challenge, here 
we restrict our analysis to such low number densities p 
of the colloids that the mean distance p~ x l d between the 
colloids is large compared with the range of the CCFs, 
i.e., p~ x l d > 2R + £. In this limit the approximation of 
pairwisc additive CCFs is expected to be valid. 



For the pair potential V c 
scaling theory predicts 



id) 



k B TU^ d) of the CCF 



C/ C (d) (£>; R, t, h b ) = (6, A, A) , (21) 

where 0^ is a universal scaling function and the scal- 
ing variables are O = sgn(t) D/£, A = D/R, and 
A = sgn(h b )D/^ h \ Within MFT for the Hamilto- 



-1/3 



Uc^ depends on the sign of 

A because the surface fields h s imposed by the colloids 
break the bulk symmetry w.r.t. hf, — > — h^. For the fol- 
lowing study it is useful to introduce another scaling vari- 
able E = A/ |0|, which depends solely on the properties 
of the solvent but not on the surface-to-surface distance 
D between the colloids. Therefore, Uc is given by an- 
other, also universal scaling function (compare Eq. ([3])) 



nian given in Eq. (|2U1) one hasi^ £ (t > 0) = |r| and 

£(/>) = 3-1/2 



Ui d) (D; R, t, h b ) = A-V d) (0, A, E) 



(22) 



with the scaling variable E = sgn (hb)£,/^ h '. 

In d = 3 the scaling function 6*( d=3 ) is not known for 
the full range of the scaling variable A and there are no 
results concerning its dependence on E. However, if the 
colloid radius is large compared to £ only those surface- 
to-surface distances D matter which are small compared 
with R. This implies the validity of the Derjaguin ap- 
proximation which allows us to express 9^ in terms of 
the universal scaling function for the film geometry. The 
latter is known from MC simulation a 23 ' 24 in d = 3 and for 
hb = 0. In the opposite, so-called protein limit D ;§> R 
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additional knowledge about the CCFs is available! 17 ' 18 
For the configuration of a colloid near a planar wall it 
has been found^i^i^ that the Derjaguin approximation is 
valid up to A = D/R < 0.3. Concerning a discussion of 
the Derjaguin approximation for the sphere-sphere geom- 
etry see Ref. [l7| and in the case of non-spherical objects 
near a wall see Ref. [l9|. For the configuration of two 
spheres as studied here one ha a 8 ' 17 ' 30 

/oo 
(x- 2 -x- 3 )^(xQ,X) dx, 

(23) 

where t?| (y, £) is the universal scaling func- 

(s) 

tion of the CCF Jq\\ per kgTc and per area 
A for a slab of the thickness L: fc,\\/A = 

k B T^L- d d { ^ (y = sgn (t) L/C S). We point out 
that within the Derjaguin approximation the scaling 
function 6^ d=3 ) does not depend on A, which therefore 
enters into Eq. ([2"2"]1 only as a prefactor. 

There are experimental indication o 63 ' 64 and theoreti- 
cal evidence^ 8 , for a pronounced dependence of CCFs on 
the composition of the binary liquid mixture acting as a 
solvent. As discussed in Appendix |A1 this translates into 
the dependence on the scaling variable S. In order to be 
able to capture this dependence to a certain extent on 
the basis of presently available theoretical knowledge we 
propose the following approximation for (y, £): 

<' fe - £) -< >fos ° ^pI|s!o) - (24) 

This approximation offers three advantages: (i) For d —> 
4, i.e., for MFT, the rhs of Eq. reduces to the cor- 
rect expression for the full ranges of all scaling vari- 
ables, (ii) For h b — > the rhs of Eq. (|24|) reduces ex- 
actly to i?y (y, £ = 0) for all d. In a certain sense the 
MFT approximation is concentrated in the dependence 
on hi,, (iii) The MFT treatment of the dependence on 



hb does not suffer from not knowing the amplitude of 
i?|l <i_4 '' within this approximation; for d < 4 on the rhs 
of Eq. ([24j) this amplitude drops out. Within the MFT 
expressions in Eq. (|24[) the scaling variables O (entering 
via Eq. (|23|) ) and £ are taken to involve the critical bulk 
exponents in spatial dimension d so that the approxima- 
tion concerns only the shape of the scaling function itself 
which typically depends on d only mildly (see, e.g., the 
comparison of MFT results with corresponding results 
obtained by MC simulations for d = 3 or with exact 
results for d = 2 in Refs. [23| and l62h. Here we have cal- 
culated S* =4) also for £ ^ 0. For ^ =3) (y, E = 0) we 

use the MC simulation results.— The scaling functions 
resulting from Eq. and from the local functional 

metho d 21 ' 22 are comparable^ and in qualitative agree- 
ment with corresponding curves provided in Ref. 1271 

It is most suitable to determine the CCF from §mft 
via the so-called stress-tensor The CCF per area in a 
slab, which is confined along the ^-direction, is given by 
the (z, z) component of the thermally averaged stress ten- 
sor, f^ FT) /A = k B T^ {T z ,z [i^MFt] - T z ,z [<t>b,MFT]), 

with 

Tz, z W\ = ^ {4>'f ~ ^4? - ||0 4 + h b cf), (25) 

where the OP <\> — <f> (z = zq) and its derivative 4>' = 
(d(j)/dz) z=Za arc evaluated at an arbitrary point —L/2 < 
zo < L/2 within the slab and 0t is the bulk OP. (The 
so-called "improvement" term of the canonical stress- 
tensor— can be neglected because it does not contribute 
to the CCF.) 

III. RESULTS 

A. Range of parameters 

In our study we use the following pair potential U (see 
Eqs. Q and ([23|)): 



U{r = D + 2R) 



| oo 
\s 



o x < 

{aexp(-x) + (l/aO^J (*/££)} x > 0, 



(26) 



where s = kR, a = A/s, ( = sgn (t) k£, and x = kD = 
nr — 2s. (The amplitude a should not be confused with 
the acronym for the preferentially adsorbed phase.) This 
parametrization has the advantage, that the shape of U 
is determined by a, C, and S, while s tunes the overall 
strength of the potential without affecting its shape. The 
ratio of the competing length scales of repulsion and of 
the CCF are measured by Cj which is typically varied 
experimentally; a is usually kept constant and provides 



a measure of the repulsion, while S (Eq. (|22|) ) depends 
solely on the thermodynamic state of the solvent. 

In the following we discuss the ranges of the values 
of the parameters entering into the effective potential 
(Eq. ([2^]) ) and of the scaling variable E which corre- 
spond to possible experimental realizations. The radius 
R of colloidal particles typically varies between 0.1/xm 
and l/i7«. In the experiments reported in the present con- 
text so far colloidal suspensions had been stabilized by 
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electrostatic repulsion, with the value of the strength pa- 
rameter A ranging over several orders of magnitudes, i.e., 
from A ~ 10 2 up to 10 5 . . . 10 6 (see, e.g., Refs. 0, B and 
[47h . The value of can be tuned by salting the so- 
lution. Due to screening effects, an increased amount 
of ions in the solution leads to a decrease of Al- 
though the coupling between the charge density and 
OP fluctuations is not yet fully understood, there is 
experimental^ 7 - and theoretica l 6 ^ 69 evidence that criti- 
cal adsorption and CCFs can be altered significantly by 
adding ions to the binary liquid mixture. Such subtle 
mechanisms are not taken into account within the effec- 
tive potential discussed here (Eq. ([26)0 . Therefore it is 
applicable only for not too small values of the screening 
length, i.e., for k" 1 > lOnm. For binary liquid mixtures 
the correlation length amplitude £0 is of the order of few 
Angstrom. The relevant experiments have been carried 

(s) 

out at room temperature due to T c ~ 300.fr. In those 
experiments deviations from the critical temperature as 
small as T — ~ 10mK have been resolved^ which 
corresponds to a correlation length £ of a couple of tens 
of nm. 

It is more difficult to assess the experimentally relevant 
range of the scaling variable S, which is a function of the 
bulk ordering field hj,. Often the amplitude Q of the 
correlation length ^ h > (hf,) is not known. However, one 
may use the equation of state (Eq. (|A4[) ) which relates 
S to the scaling variable X associated with the OP 4> 
(see the text after Eq. (|A3j) ) . In terms of the parameters 
of the potential given in Eq. (|2"6"|) . for t > one has 
X = moC~ 1//,y with (see Appendix [A"|) 

m = Sff n(0)(C o + ) 1/l/ IS/01 1/(3 (27) 

and £q~ = re£g~. For example, in binary liquid mixtures 
the order parameter </> is proportional to the deviation of 
the concentration c a of the component a from its criti- 
cal value ci'Jc, 4> = A(c a — c^l) (note that B is propor- 
tional to A), which can be easily controlled by chang- 
ing the mass or the volume fraction of one of the com- 
ponents of the mixture. The experiments reported in 
Refs. and HI provide indications concerning the size of 
the critical region in the thermodynamic direction or- 
thogonal to the temperature axis for the binary liquid 
mixture of water and lutidine near the consolute point 
of its phase segregation. These measurements revealed 
the occurrence of CCFs within the range of the luti- 
dine mass fraction deviating from its critical value 
w l.c up to \u>l —ojl.c\ — 0.04. From the experimental 
data in Refs. [70I - I72I one findsZ 3 - for the water-lutidine 
mixture B u « 1.0. The index u) refers to the spe- 
cific choice of the order parameter, i.e., ^ (f 0") = 
ujl — k>L,c = B u \t\ . Fitting the experimentally deter- 
mined coexistence curv o 47 ' 63 yields a somewhat smaller 
value £> w = 0.765 which we adopt in the following. There- 
fore the difference w l,c| — 0.04 corresponds to 
\B/4>\ — 19, which for, e.g., (q = 0.02 translates into 
I mo I ~ 17 (Eq. (|27p). where we have used the criti- 



cal exponents of the three-dimensional Ising universality 
class4£ v(d = 3)= 0.6301(4) and /? (d = 3) = 0.3256(3). 
(The critical composition corresponds to mo = ±00.) Ac- 
cordingly, for the temperature differences accessible in 

these experiments, i.e., for (t - Tj s) ) /T c (s) w 3 x 10~ 5 , 

one has for the scaling variable \X\ « 0.26, corresponding 
to |E| w 6.6. 

As discussed before, the effective pair potential given in 
Eq. (|26|) is applicable only for sufficiently large distances 
D > k^ 1 because it takes only the interactions of the 
double-layers into account and neglects possible short- 
ranged contributions to effective van-der-Waals interac- 
tions. Furthermore, the critical Casimir potential takes 
its universal form (Eq. (|2"2"|) ) only in the scaling limit, 
i.e., for distances D which are sufficiently large compared 
with the correlation length amplitude £q~ ~ 0.25nm. 
Analogously also £ and £( h > must be sufficiently large 
compared with microscopic scales. Later on, in or- 
der to circumvent the unphysical divergence ~ x~ l in 
U (D — > 0) (see Eq. (|26|) ) for small distances D we shall 
consider, as far as necessary, also a linear extrapolation: 

U(D< D ) = Acxp(- K D) + U cfi + (D - D )U' cfi , (28) 

where D ~ £ , U c = U C {D = Dq), and U' c0 = 
(dU c /dD) D=D . 



B. Thermodynamics 

In this section, we consider pair potentials which are 
repulsive at short distances (i.e., a is sufficiently large) 
so that the suspension is stable. 



1. General discussion 

In this Subsection we consider the thermodynamics of 
actually ternary colloidal suspensions with binary sol- 
vents such as water-lutidine mixtures which exhibit a 
closed-loop two-phase region of demixed phases (each 
being rich in one of the two species). We focus on 
that region of this miscibility gap which is close to the 
lower critical point. For fixed pressure, their thermo- 
dynamic states can be characterized by the tempera- 
ture T and the concentration c a of one of the species 

with the critical point (^Tc^ , c^cj and the liquid-liquid 

phase coeaistence curve Tex (c a ) in the absence of col- 
loids. Upon adding colloidal particles to such a solvent, 
for fixed pressure the thermodynamic space of the system 
becomes three-dimensional spanned by T, c a , and, e.g., 
by the colloidal number density p (see Fig.[TJ one can also 
choose, instead, the fugacity of the colloids). Accord- 
ingly, the closed-loop phase coexistence curve Tcx{c a ) 
becomes a two-dimensional, tubelike manifold T cx (c a , p) 
with T cx (c a ,p = 0) = Tcx(c a ). It contains a line C c 
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FIG. 1. Sketch of the phase diagram for colloids immersed in a binary liquid mixture, the latter exhibiting a closed-loop 
miscibility gap. We focus on the region around the lower critical point 13 (Tc , ci^c, p = 0) (see the lavender phase separation 
curve Tex 1 (c a ) in the (T,c a ,p — 0) plane). For fixed pressure corresponding to a liquid state of the system, the thermodynamic 
(td) space of the pure solvent is two-dimensional; here we consider as variables the temperature T and the concentration 
c a . Upon adding colloids at fixed pressure the td space becomes three-dimensional; we use the colloidal number density p as 
additional td variable. The miscibility gap of the pure solvent in the (T, c a , p = 0)-plane extends to a two-phase region in the 
three-dimensional td space and is bounded by a two-dimensional manifold T cx (c a ,p) of coerristing states (not shown). The 
tube-like shape of T cx (c a ,p) is not straight but bent and twisted. Its actual form is expected to depend sensitively on all 
interactions (i.e., the solvent-solvent, the solvent-colloid, and the colloid-colloid interactions). Each state on T cx (c a ,p) coexists 
with another one, both being connected by a horizontal and straight, so-called tie line (green dash-dotted line). There is a line 
C c (black curve) of critical points (some of which are shown as black squares) embedded in T cx (c a ,p) which is the extension 

of the critical point ^Tj s ', ci°', p = OJ. C c is given by the states for which the tie line has zero length. The projections of C c 

(with selected critical points ■) onto the planes (p, c a ) and (T, c a ) (with + as the projection of ■) are indicated as red dashed 
curves. For p — > the bending of the curve C c can be inferred from scaling arguments (Eq. (|29[) ~). For colloidal suspensions 
interacting via an effective potential which consists of a soft repulsion and the critical Casimir attraction, the results of the 
effective approach suggest that, for intermediate values of p, C c bends up again (see Subsect. HUB 4[) . This bending is due to 
the specific properties of the critical Casimir forces and is indicated in (a). Effective colloidal models render coexistence curves 
Tcx* {p | c a ) which are explicit functions of p only and depend parametrically on the overall concentration c a (see (a) for three 
examples; the three vertical arrows indicate thermodynamic paths describing the approach of the corresponding critical point ■ 
upon raising temperature). For suitable effective models and within a certain region of the td space (see main text), the critical 

points ( x), (Tc"*^ (c ) , p c (ca) j , of the coexistence curves Tc X **\p \ c a ) are expected to approximate the projection of C c onto 

the (T, p)-plane. Within the effective approach a unique value c a is taken throughout the whole system. In contrast, all three 
panels show that in general for T — const the coexisting phases (i.e., the points connected by a tie-line) differ both in p and c a - 
Thus the effective approach has a limited applicability for determining the phase diagram. Experimentally, upon increasing 
temperature one is able to determine a (in general nonplanar) coexistence curve (black line in (c)) in the three-dimensional td 
space. In (c), for a selected critical value f) c = (1 — c a , c (p) , p) (i.e., f) c is a point on the red dashed line in the (p, c a )-plane) such 
a curve (black line), its projection onto the (T, p)-plane (mustard line), as well as a selected tie line and its three projections 
are shown. Note that the tie line is not parallel to the p-axis. This means that the coexisting phases differ with respect to 
c a , in contrast to the aforementioned assumption of an unique value of c a for the effective one-component description. Thus 
the mustard curve in (c) will in general differ from the corresponding curve T^x '{p \ c a ) in (a), even if the associated critical 
points are the same. 
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of critical points (T c (p) , c a , c (p)) which is the exten- 
sion of the critical point of the solvent in the absence 



of colloids, i.e., (T c (p - 
For fixed temperature 

(2) 



0),c QiC O = 0)) = (T, 



(T c (s) ,c^) 



c (1) 

<-a , P 



(2) 



T = const, the set of pairs 
of coexisting states given by 

T = T cx (c a , p) forms a curve C cx (T) in the (c a , p) plane, 
the shape of which depends on the considered value of T 
(see Figs.QJa) and (b)). In the phase diagram, two coex- 
isting states are connected by a straight, so-called tie line 
(see Fig. []}. (If one chooses the fugacity of the colloids 
instead of p, one also obtains a tubelike manifold T cx of 
phase coexistence. However, in this case the horizontal 
tie lines lie in the plane of constant fugacity, i.e., parallel 
to the (T, c a )-plane, because the coexisting phases share 
a common fugacity of the colloids.) 

As stated above, the two-phase loop of the pure sol- 
vent (i.e., for p = 0, bounded by Tcx(c a )) extends into 
the three-dimensional thermodynamic space of the ac- 
tual colloidal suspension (p ^ 0). Due to the presence 
of additional interactions and degrees of freedom one ex- 
pects that the shape of this two-phase region (bounded 
by T cx (c a ,p)) is not a straight but a distorted tube. The 
actual shape of T cx (c a , p) is expected to depend sensi- 
tively on all interactions present in the ternary mixture, 
i.e., among the colloidal particles, between the colloidal 
and the solvent particles, and among the solvent parti- 
cles. The relevance of the solvent-solvent interaction for 
the effective potential and, accordingly, for the phase be- 
havior of the effective colloidal system has been demon- 
strated recently by MC studies in which various kinds of 
model solvents have been used^ It is reasonable to ex- 
pect that this relevance transfers also to the phase behav- 
ior of the full multi-component system. Such distortions 
of the phase diagram relative to that of the underlying 
binary mixture do occur for ternary mixtures of molecu- 
lar fluids. For example, in Ref. [T5| experimental studies 
of molecular ternary mixtures containing various kinds of 
lutidines are reported. These studies show that the up- 
per and lower critical temperature for a closed-loop phase 
diagram can be tuned by varying only the concentration 
of the third component and that the two-phase loop can 
even disappear upon adding a third component. Similar 
experimental results are reported in Ref. l76l . Such com- 
plex phase diagrams can also be generated by adding 
colloidal particles to the binary solvent, as can be in- 
ferred from corresponding experimental studies^ - — In 
contrast to the molecular ternary mixtures, for the latter 
kind of ternary mixtures a decrease of the lower critical 
temperature upon adding colloids as a third component 
is reported. 

Theoretical studies of bona fide ternary mixtures have 
so far been concerned with, e.g., lattice gas models 7 - 7 - and 
(additive or non-additive) mixtures of hard spheres, nee- 
dles, and polymers . 78 ' 79 In these studies the constituents 
are of comparable size, i.e., their size ratios are less than 
ten. The peculiarity of the kind of mixtures considered 
here lies in the fact that the sizes of their constituents 



differ by a few orders of magnitude. This property dis- 
tinguishes them significantly from mixtures of molecu- 
lar fluids. In contrast to molecular ternary mixtures, in 
colloidal suspensions the colloidal particles influence the 
other two components not only by direct interactions but 
also via strong entropic effects. This is the case because 
their surfaces act as confinements to fluctuations of the 
concentration of the solvent and they also generate an 
excluded volume for the solvent particles. The impor- 
tance of considering the colloidal suspension as a truely 
ternary mixture has been already pointed out in Ref. l35l . 



2. Scaling of the critical point shift 

For dilute suspensions, i.e., for p — > 0, the shape of 
the line C c of the critical points can be estimated by re- 
sorting to phcnomcnological scaling arguments similar to 
the ones given by Fisher and Nakanishi 8 ^ for a critical 
medium confined between two parallel plates separated 
by a distance L. For a dilute suspension, the mean dis- 
tance p~ x l d between colloidal particles plays a role anal- 
ogous to L. Close to the critical point of the solvent, 
due to £ ~ \t\ " one can identify the two relevant scal- 
ing variables oc \t\~ hb and cx \t\ u p~ x / d (for 
the simplicity of the argument 8 -^ here we do not consider 
the influence of the scaling variable \t\ v R) and propose 
the scaling property f (T, Ap,p) ~ \t\ ~ a K (w^- 1 ' ,vj( 2 ') 
of the free energy density (where the difference A/i ~ hb 
of the chemical potentials of the two components of the 
solvent acts as a symmetry breaking bulk field). The 
critical points are given by singularities in K occurring 

at certain points \ Wc \ wi 2 ^ . This implies that the crit- 
ical point (T c (p) , hb lC (p)) shifts according to 



T c (p) - T« ~ p x l {vd 
so that 

Ca.c (p) 



and h b ,c (p) ~ p A/( " d) 

(29a) 



A/(dSv) 



(29b) 



with (1/ M)) d=3 ^ 0.53 and (A/ (dSu)) d=3 ~ 0.17. This 
states that in the presence of colloids the critical point 
occurs when the bulk correlation lengths £ and t;^ of 
the solvent become comparable with the mean distance 
p-i/d between the colloids. 



3. Effective one-component approach 

Integrating out the degrees of freedom associated with 
the smallest components of the solution (here two) pro- 
vides a manageable effective description of colloidal sus- 
pensions. This kind of approach is commonly used, for in- 
stance, recently for studying large particles immersed in 
various kinds of model solvents^ or in order to describe a 
binary mixture of colloids immersed in a phase separating 
solvent* 8 ^ However, this effective description has only a 
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limited range of applicability for investigating the phase 
behavior of colloidal suspensions. For example, it fails in 
cases in which the influence of a set of colloids on the sol- 
vent cannot be neglected or the (pure) solvent undergoes 
a phase seperation on its own, as considered here. In the 
effective approach, the concentration c a of one of the sol- 
vent species, averaged over the whole sample, enters as a 
unique parameter into the effective interaction potential 
between the large particles (see the dependence on S in 
Eqs. ([212]) and (|2U]0 . However, generically (e.g., due to 
the adsorption preferences of the colloids) the COnCentra- 
fi) 

tions c a arc different in the two coexisting phases i = 1,2 
(see Fig. [I]) and thus, for a proper description, one would 
have to allow this parameter, and hence the effective po- 
tential, to vary in space. Experiments have revealed^r— 
that for suspensions very dilute in colloids with a phase 
separated solvent, basically all colloidal particles are pop- 
ulating the phase rich in the component preferred by the 
colloids (with concentration c a )■ In this case the ac- 
tual effective attractive interaction among the particles 
will be weaker than implied by the effective potential 
in which only the overall concentration c a < enters 
as a parameter. (We recall that the CCFs depend non- 

(s) 

monotonously on c a and are strongest for c a < c a . c . Thus 

for S a 2 ^ < c a < c a the CCFs differ for each concentration 
and are - for the typical situation considered here - in gen- 
eral weakest for c a ^ ; c a 2 ^ is the concentration in the other 
coexisting phase.) Within the effective approach for the 
colloids, one obtains, e.g., by means of DFT, coexistence 
curves Tcx^\p \ c a ) which depend parametrically on the 
solvent composition c a « c a % (or equivalently, by making 
use of Eq. (|A"4|) . on £ as in Eq. ^Q). Since the effec- 
tive potential corresponding to this specific value of c a 
is used both for the one-phase region as well as for both 
phases within the two-phase region, this implies that the 
tacitly assumed corresponding physical situation is such 
that the composition c a is fixed throughout the system as 
an external constraint. In particular the two coexisting 
phases do not differ in their values of the concentration 
of the solvent particles but only in their colloidal densi- 
ties p^ 1 ' and p( 2 >. Thus within the effective approach for 
determining the phase behavior one of the essential fea- 
tures, i.e., the tendency of the solvent to phase separate, 
is suppressed. Accordingly, it is not possible to construct 
the full coexistence manifold T cx (c a , p) on the basis of 

the curves Tex (p I c a ) alone. Rather, the effective ap- 
proach is adequate as long as to a large extent the phase 
segregation involves only the colloidal degree of freedom, 
i.e., the values of p differ in the two phases, but the values 
of c a are nearly the same. Therefore the approximation 

T cx (c a , p) ~ 7«c (p I c a) is valid in a region of the ther- 
modynamic space in which the actual tie lines happen to 
be almost orthogonal to the c a -axis (see Fig.Q}. For the 
phase diagram in terms of the variables T, c a , and the 
fugacity of the colloids this latter condition implies that 
the tie lines have to be sufficiently short. 



We expect the effective one-component approach to 
work well for temperatures corresponding to the one- 
phase region of the pure solvent and for an intermedi- 
ate range of values of p. On one hand p should be large 
enough so that the competition between the configura- 
tional entropy and the potential energy due to the ef- 
fective forces can drive a phase separation. On the other 
hand p has to be small enough so that the approximation 
of using an effective pair potential between the colloids is 
valid and the influence of the colloids on the phase behav- 
ior of the solvent is subdominant. Given such values of p 
and T, the reduced second virial coefficient (Eq. Q14p) 
is an appropriate measure for the strength of the attrac- 
tion and a useful indicator of the occurrence of a phase 
separation into a colloidal-rich ( "liquid" ) and a colloidal- 
poor ("gas") phase. According to the discussion above, 
for values of the thermodynamic variables (i.e., the val- 
ues of T and p, and for the prescribed value of c a ) which 
are approximately the same as in the one-component 
approach, the ternary mixture is expected to exhibit a 
phase separation. According to the VLNF conjecture 
the corresponding critical point should occur when i?| 
reaches the critical value £?2. c — -1.212 (Eq. dU])) of 
Baxter's model. While the VLNF conjecture provides an 
empirical estimate for the values of parameters for which 
there is a critical point, within DFT these parameters as 
well as the shape of the phase coexistence curve can be 
calculated. Within the effective approach only its depen- 
dence on the colloidal density number p for globally fixed 
values of c a can be determined. 



4. Phase diagrams 

Generically, in experiments the solvent composition 
and hence the associated quantity m (Eq. (|27|l ) is fixed 
and £ ~ \t\ v is varied. In Fig.[2ja) for a solvent exhibit- 
ing a lower critical point and for the parameter choices 
aT/T c (s) = 100, s = 10, m = -100, -20, and -6, and 
Co = 0.01, 0.05, and 0.1 we present the coexistence curves 
Tcx^\p | c a ) and the spinodals of the colloidal "liquid" - 
"gas" phase transition as function of the colloid packing 
fraction n = (4ir/3)R 3 p as obtained by the DFT pre- 
sented in Eq. (fT8)l . (Note that in accordance with Eqs. ([2]) 
and (|26|) in the product aT the explicit dependence on 
T drops out.) For the corresponding free energy J-rpa 
given in Eq. (fT"9| the effective colloidal system phase sep- 
arates if U a0 < (u a ,oj = —21.3. This condition is 

satisfied provided the attractive part of the interaction 
potential is sufficiently strong. We recall, that for the 
effective potential considered here attraction occurs for 
T/r c (s) = 1 ± (C/Co)~ 1/!/ -> 1 (see Eq. ([26])). For a sol- 
vent exhibiting a lower (upper) critical point the lower 
(upper) sign holds in the one-phase region of the sol- 
vent, within which the effective approach is applicable. 
In this temperature limit the variation of a (Eqs. @ and 
(|26p ) with T is subdominant and thus the dependence of 
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FIG. 2. (a) Phase coexistence curves T<£ (p \ c a ) (full lines), spinodals (mean-field divergence of \T> dotted lines), and the 
critical points Tc (crosses) of an effective, one-component system of large particles (Eq. (J2H|)) as obtained by DFT according 
to Eq. (|18[) . The curves correspond to a solvent with a lower critical temperature T"i s ' and to various fixed solvent compositions 
c a {mo = —100 (red), —20 (green), —6 (blue)). The other parameters are taken as £o = 0.01 (□), 0.05 (o), and 0.1 (A), s = 10, 
and aT/Tc = 100. Close to the phase separation, the dominant temperature dependence is that of the critical Casimir forces 
(CCFs) encoded in £ ( t = 1 — T/tJ s ' J = sgn (t) s£ (t) = sgn (t) Co Therefore, in terms of C (see panel (b)) the curves in 



(a) for different values of Co fall de facto on top of each other, independent of whether the solvent has an upper or a lower critical 
point. Within the integral equation approach only the spinodals have been determined (PY: dashed lines in (b), see the main 

text). According to Eq. (|27[) . [c a — c a s ,Vj ~ mo — ±oo corresponds to the critical composition c a B } c , For solvent compositions 
which are somewhat poor in the component preferred by the colloids, i.e., for intermediate negative values of mo — —20, the 
CCFs are strongly attractive. Therefore, for them small correlation lengths suffice to bring about phase separation; accordingly 
the binodals occur at small values of C ~ \t\~" ■ Here we consider only thermodynamic states of the solvent which are in the 
one-phase region, i.e., t > 0. 



Frpa on T /Tc and Co reduces to a dependence on C 
only. Accordingly, as can be inferred from the compari- 
son of Figs. r2fa) anddfb), in terms of C for each value of 
mo the coexistence curves Tcx^\p \ c a ) for different val- 
ues of Co fall de facto on top of each other. The difference 
between these curves is of the order of (Co/C) 1 ^ (because 
T/T c (s) = 1 ± {(/( y 1,u ), which for the three values of 
Co used in Figgis about the thickness of the lines shown 
in Fig. [UJb). In terms of this presentation it does not 
matter whether the solvent exhibits a lower or an upper 
critical point. Although spinodals are mean- field arti- 
facts, we present them nonetheless because they provide 
some indication about the location of the binodal which 
encloses the former. Furthermore the spinodals carry the 
advantage that the isothermal compressibility \t is the 
property of only one thermodynamic state. Therefore, 
in contrast to the calculation of the binodal (which de- 
pends on two coexisting phases), the calculation of the 



spinodal does not suffer from the non-uniqueness of the 
effective potential in the case of phase coexistence. As 
discussed in Subsect. Ill A 21 based on formally exact re- 
lations the phase behavior can in principle be calculated 
from the correlation functions obtained within the inte- 
gral equation approach (IE A). Within the so-called com- 
pressibility route (Eq. ©) the spinodals, i.e., the loci 
of the mean- field divergence of \t, are directly accessi- 
ble. On the other hand, the binodals, i.e., the loci of 
two thermodynamic states which at the same tempera- 
ture have different packing fractions but the same pres- 
sure, are directly accessible via the so-called virial route 
(Eq. (|4|)). We refrain from calculating the binodal (the 
spinodal) via the compressibility route (virial route) , be- 
cause it would require thermodynamic integration, which 
we want to avoid due to the subtlenesses described in 
Subsect. Ill A 21 The spinodals calculated by the com- 
pressibility route are shown as dashed curves in Fig.[^b). 
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For the region of the thermodynamic space where the 
IEA renders solutions for the correlation functions, we 
failed to find binodals, i.e., along the pressure isotherms 
p (77; T = const) as calculated via the virial route there 
are no two states with 771 ^ 772 and p (771; T) = p (772; T). 
These observations can be explained by the thermody- 
namic inconsistency of this IEA. Due to the approximate 
bridge function, the binodals as obtained by two different 
routes need not to coincide. The same observations are 
found within the HNC approximation within which the 
spinodals are shifted w.r.t. to the PY results to slightly 
smaller values of (. Although within the IEA the bin- 
odals could not be determined, at least the loci of the 
obtained spinodals are similar to the ones obtained from 
DFT. 

Since the CCFs are strongest for slightly off-critical 

compositions 777^ 1 ~ ^c a — c£}cj < the binodals (and 

spinodals) are shifted to smaller values of C upon decreas- 
ing |m | down to m ~ —20. With a further decrease 
of I mo I the system moves too far away from the critical 
point of the solvent so that the CCF weakens and the 
spinodals shift again to larger values of £. 

The critical value 77 c of the packing fraction is rather 
small, i.e., ?7 C w 0.07 (see Fig. [2]), because the effective 
hard sphere diameter a which results from the soft repul- 
sion U rep is larger than 2R. In terms of n a = (a/(2R)) 3 rj 
(defined after Eq. (|19p ) the critical value assumes its 

RP A- value rji R c PA ^ = 0.129. Furthermore, the binodals 
shown in Fig. [2] are rather flat compared with, e.g., the 
ones for hard spheres interacting via a short-ranged, at- 
tractive temperature independent potential (described 
also by Eq. (fl"9j) ). In the present system, the devia- 
tion from the critical temperature T c which leads to a 
range of 77 for the coexisting phases as large as shown 
in Fig. [21 is about l%o, whereas for a system of hard 
spheres with an attraction the corresponding tempera- 
ture deviation is a few percent. For smaller values of (0 
the binodals are flatter and the differences of the critical 
temperatures for different values of mo are smaller (see 

Fig. Ufa)). T he requ irement 77- 1 / 3 > (4tt/3) 1/3 (2 + (/s) 
(see Subsect. IIIB2j) concerning the validity of the pair- 
wise approximation for the CCFs is fulfilled for the whole 
range of values for 77 and £ shown in Fig.[2j e.g., for s = 10 
and C = 4 the above condition is 77 < 0.3. 

The sketch of the phase diagram in Fig. Q] and the cal- 
culated coexistence curves in Fig. [2] correspond to using 
the background potential with fixed parameters s and 
a. In order to investigate the dependence of the locus 
of the phase separation on properties of the background 
potential, in the following we shall discuss how the criti- 
cal temperature T c (e//) = T c (e//) (s, a, s, Co, T^ s) , ±) of 

the effective colloidal system varies as function of its 
arguments (± is taken to be — and -I- for a lower and 
an upper critical point, respectively). According to the 
discussion above, to this end it is sufficient to deter- 
mine Ce 6 ^ (S,a, s). In Fig. [3] for fixed values of s and 
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FIG. 3. Attraction due to CCFs can induce a "liquid"- 
"gas" phase separation of the colloids. Within the effective 
approach described by Eq. (I26|l the dominant temperature 
dependence is that of the Casimir scaling function encoded 
in £ — K^(t) ~ \t\ ~ v . The full curves correspond to DFT 
(Eq. (JTHJ) and show for which values of E and a correspond- 
ing to the solvent composition and the repulsive strength of 
the background potential, respectively, a given temperature 
(in terms of Q is the critical one Tc"^^ or Ci e ^' . The contour 
lines a (S; s; Cc"^') depend parametrically on s = kR where 

k" 1 is the range of the background repulsion and R is the 
radius of the colloids. Systems, with their parameters being 
located in the area below the contour line (^c e ^\ are phase 
separated at this temperature. This two-phase region expands 
to larger values of a with increasing For a fixed value Q e , 
a (E; s; C,i e!! ^ is largest for compositions slightly poor in the 

component preferred by colloids, i.e., E < 0. In addition, 
the corresponding contour lines based on the VLNF criterion 

B2 (c = Cc" 5 ^^', E, a; s) = 13% <c for the reduced second virial 

coefficient (Eq. (|14p 1 are plotted as dashed lines. It is satisfac- 
tory to see that both approaches provide comparable results. 



for specific values of Ce 6 ^ we show the contour lines 
a = a (S; s; C^^^ ■ That is, for a colloidal system with 
a background potential characterized by certain values 
of a and s, from the plot in Fig. [3] one can read off 
the value of the difference of the chemical potentials 
of the two components of the solvent (which is related 

to S) for which Tc & ^^ has a prescribed value. The 
full curves are calculated within the density functional 
approach (Eq. (fT5|) ) and are compared with the ones 
corresponding to the simple prediction (Eqs. (| 13[) and 
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E, a; s 
Wf) 



B*2, c 



-1.212, which 



also yields a relation C,c CJ J> = Cc eJ ' M (S,a;s). Both ap- 
proaches differ only slightly. For given values of Cc 6 ^' 
and E the VLNF conjecture predicts a smaller critical 
value of a than obtained from DFT. The largest devia- 
tions occur for slightly negative values of E, for which the 
CCFs are more attractive. For all values of the param- 
eters shown in Fig. [3] (i.e., a > 30), at small distances 
D the corresponding potentials exhibit repulsive barri- 
ers > 30fcsT. For smaller values of a, i.e., a < 15, this 
potential barrier disappears and the colloidal particles 
may form aggregates. This is discussed in the following 
paper^ see Figs. 1 and 2 therein. 



The contour lines a = 



E;,;d e//) 



attain their 



largest values for slightly negative values of E, i.e., for 
E ~ —2. This again reflects the asymmetry of the 
strength of the CCFs with respect to the bulk order- 
ing field hb in the presence of symmetry breaking surface 
fields. In Fig. [3l for all parameter pairs (a, E) below each 
contour line, characterized by (d ; the corresponding 
system at the temperature belonging to the value ci e ^ is 
phase separated. This two-phase region widens with in- 
creasing values of £ and it significantly expands to larger 
values of a, which is in line with the fact that the CCFs 

(s) 

become stronger upon approaching T c . For s = 10, 
upon increasing C^^^ from 3 to 5 the largest critical 
value of a increases by two orders of magnitude (com- 
pare the corresponding blue and green curves in Fig. [3]). 
The two-phase region can be markedly increased upon 
increasing s (compare the blue and gray curves in Fig. [3] 
corresponding to Q = 3 with s = 10 and s = 50, re- 
spectively) because for the phase separation the strength 
of the attraction is important; due to U ~ s (Eq. (|26|) ) 
the latter can be easily tuned by varying s. 

For the experimental studies reported in Refs. [lH and 
l33l colloidal mixtures consisting of silica spheres, water, 
and 2-butoxyethanol or lutidine have been used; with 
a radius R w llnm the colloidal particles have been 
rather small. Thus, a priori, the effective approach, de- 
veloped in Subsect. Ill Al and discussed in Subsect. MI B 31 
is not expected to apply. Nonetheless, one can try to 
do so and estimate the critical temperature for solu- 
tions containing such small particles by using the gen- 
eral scaling arguments presented in Subsect. MI B 21 and 
in Ref. At the critical temperature, the dimension- 



,(3) 



t\ 



R associated with the 
with 



less scaling variable 

(3) 

radius of the colloid takes on a certain value uu c 

(3) / (: 
Wc '/wfr 

(3) 



(T c (s) - T c {pj) /T c (s) | R = \t c \"R (where 







is a non- universal amplitude). 



tive approach we can calculate w c 3 ^ /w^ 



Using our effec- 
within mean- 
field theory for parameters corresponding to the solvents 
used in those experiments, and demand^ that this ra- 
tio applies approximately to the experiments reported 
in Refs. [H and For the solvents used, the relevant 



parameters are £0 ~ 0.25nm and £0 ~ 0.03. Choos- 
ing s = 10 (see Subsect. IIII Al and Fig. Ufa)) renders 
R = Cos/Co ~ 85nm. This value of R is sufficiently 
large so that the aforementioned effective approach can 
be considered to be reliable. From Fig. [2ja) we can read 
off the corresponding shift \t c \ ~ 0.001 which leads to 

w c 3 ^ /w^ ~ 1.094nm. By using this latter mean-field 
estimate, for the experimental value R « llnm of the 
radius of the colloidal particles one thus obtains the shift 



l*c| 



- /„„(3) 



„(3) 



-1/v 



0.234. With the critical tem- 



perature T, 



0) 



320/v of the solvent^ 3 , this leads to 



T c (p) = (l-ic)T, 



310!^. Indeed, for such a tem- 



perature phase separation in the ternary mixture has 
been observed (although even lower critical temperatures 
T c fts 300K have been reported^^ 3 .) . As stated above 
such an agreement for even small colloids could not have 
been anticipated from the outset. 



IV. SUMMARY 

We have studied the collective behavior of monodis- 
perse colloidal suspensions with near-critical solvents. 
Colloidal particles acting as cavities set the bound- 
ary conditions for the fluctuating order parame- 
ter (OP) of the solvent and perturb the OP field 
on the length scale of the bulk correlation length 



£(«-( 



= It — T, 



/T c (s) -> 



\t\ v , which diverges 
W. 



upon approaching the critical temperature T c ; v is a 
critical bulk exponent. These modifications of the OP 
and restrictions of its fluctuation spectrum result in an 
effective force acting between the colloids, known as the 
critical Casimir force (CCF). For equal boundary condi- 
tions (BCs), i.e., equal surface properties of all colloids, 
the CCF is attractive. The CCF depends on the configu- 
ration of all colloidal particles and thus is in general non- 
additive. Wc have obtained our results for an effective 
one-component system of colloids interacting via an effec- 
tive pair potential. The regular background potential has 
been taken to be a soft repulsion acting on a length scale 
ft -1 and with strength A. We have considered suspen- 
sions with medium values of the colloidal number density 
p, for which the approximation of pairwisc additive CCFs 
is valid. The scaling function of the critical Casimir po- 
tential has been taken in accordance with the Derjaguin 
approximation (Eq. (|23[)). within which the pair potential 
between two spheres is expressed in terms of the scaling 
function of the CCF in a slab. The variation of the scaling 
function in a slab with the scaling variable E associated 
with the bulk ordering field hf, (see Eq. (J3])) in spatial 
dimension d = 3 has been approximated based on Monte 
Carlo simulation data and on mean-field theory results 
(Eq. p4]) ). On this basis we have obtained the following 
main results: 

1. In Fig. [1] we have illustrated the phase behavior of 
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a ternary mixture consisting of colloidal particles 
and a binary solvent which exhibits a miscibility 
gap. We have discussed in detail the relationship 
between this full description and the one which is 
based on an effective, one-component colloidal fluid 
(Fig. P. 

2. For certain regions of the three-dimensional ther- 
modynamic space (T,c a ,p) of a ternary mixture, 
the effective one-component model is applicable for 
obtaining the onset of phase separation. For parti- 
cles interacting via the effective pair potential given 
in Eq. (|26| we have calculated the phase coexis- 
tence curve and the spinodal (i.e., the loci where 
within mean-field theory the isothermal compress- 
ibility xt diverges) both within density functional 
theory (DFT, Eq. (fl8|)) and the integral equation 
approach. In Fig. [5] the coexistence curves for var- 

(s) — 1 

ious values of c a — c a<c ~ m are shown. The 
spinodals as obtained by using the DFT approach 
are narrower and are located at smaller values of 
£ = k£ ~ \t\ than the ones obtained from the 
integral equation approach. 

3. In Fig.[3]the dependence of the critical temperature 
of the effective one-component system (expressed in 

terms of Q ) on the parameters a = A/s, E, and 
s = kR is shown. For a given value of a, Q e is 
smallest for slightly negative values of £. The crit- 
ical temperature has been calculated within DFT 
and compared with the simple prediction B^ = B^ c 
as suggested by Vliegenthart and Lekkerkerker— 
and Noro and Frenkelj 4 - 9 is the reduced second 
virial coefficient (Eqs. (fl~3)) and ([H])) and B^ c is the 
critical value of Baxter's model. Both approaches 
yield good agreement. 

To conclude, our results show that the CCF, which can 
be easily controlled by temperature T and the strength 
and range of which can be varied by changing the bulk 
ordering field hb, can induce phase separation into a 
colloidal-poor and a colloidal-rich phase. Within the ap- 
proach of using an effective potential, we have identified 
the ranges of values for the background repulsive poten- 
tial and the values of the scaling variables associated with 
the critical solvent, for which the colloids phase separate. 
We have used the approach developed here in order to 
study also the stability of colloidal suspensions in near- 
critical binary solvents^ 

Concerning further research, more complex suspen- 
sions could be considered, for example a mixture of col- 
loids with different adsorption preferences for the sol- 
vent particles. According to Gibbs' phase rule, for 
such a four-component system (a binary mixture of 
colloids in a binary solvent) at constant pressure a 
two-dimensional manifold of critical points embedded 
in a three-dimensional manifold of coexisting states in 
the (then four-dimensional) thermodynamic space can 
be expected. Recently, such a mixture was studied 



experimentally^ In this study, concerning the effective 
interaction potential of the CCFs between the colloids 
the attractive (+, +) BCs and the repulsive (+, — ) BCs 
were realized for colloidal particles of the same kind and 
for particles of different kind, respectively. In line with 
an effective DFT, upon approaching the critical temper- 
ature of the solvent there are indications that these two 
kinds of particles phase separate. 

In order to study the phase behavior of the actual 
three- or four-component mixture for the full range of 
all thermodynamic variables T, c a , and p, an analysis 
beyond the effective approach is needed. 
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Appendix A: Off-critical mixtures 

The Hamiltonian given in Eq. (|2"0|) depends on the re- 
duced temperature t and the bulk field hb, which for a 
binary liquid mixture is proportional to the deviation 
of the difference of the chemical potentials of the two 
components a and b from its critical value, i.e., hb ~ 
{pa - Mb) - (Mo.c - Mb,c) for a fixed pressure ~ p a + Mb- 
Accordingly, for fixed hb the OP varies upon changing 
t. However, in most experimental realizations t is varied 
at fixed concentrations Cj, i 6 {a, 6}, and thus the OP 
is kept constant. The OP <p, which is not uniquely de- 
fined, is related to the ordering field hb by the equation of 
state (EOS) which for a critical system takes the scaling 
formic 

hb = Vsgn(<f>)\<t>\ s F(t\B/ct>\ 1/0 y (Al) 
or equivalcntly 

4> = sgn (h b ) B \tf G (fi 1 / W« t \h b \- 1/m ) (A2) 

with universal scaling functions F and G (Y) with 

G (Y — > 0) = \Y\~ . V and B are non-universal am- 
plitudes which depend on the definition of such 
that on the coexistence curve the bulk OP follows 
4>b (t — > 0~,hb = 0) = B\t\ . V, B, and the correlation 
length amplitudes £ ' h (defined after Eq. (|2"l"|) ) arc re- 
lated to each other by universal amplitude ratios such 
that only two of them are independent. 18 ' 40 In the lowest 
order in its argument 

X = t\B/4 ,p , (A3) 

the universal scaling function F has the functional form 
F ( X ) = l+X, which captures the crossover between the 
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critical behavior at t = and at hb = 0, respectively^ 
At the coexistence curve one has t < and hb = 0; 
therefore X = —1 in agreement with <fi = B\t\^. Along 
the critical isotherm t = one has X = and thus 
<f> = sgn(hb)\hb/'D\ 1 ' S . For our purpose of fluid sys- 
tems exposed to surfaces the sign of <f> matters and the 
appropriate scaling variable is X = sgn{4>) \t\ |S/0| 1 ^' 3 . 
In terms of the scaling variables S = sgn (hb) £,/^ h ^ = 

sgn{h b ) (e + /^' l) ) \hb\ v/m \t\~" and X, Eq. dSJ) takes 
the scaling form 



sgn(E) |S 



/3<5/V 



n5/(cJ-1) / n + 



S<5/i/ 



x sff n(X)|X|-^F±(|X 



where F± (\X\) = 1 ± \X\ and ± refers to the sign of 
t. R x , Q2, , and are universal amplitude ratios.— 
We shall use these expressions in order to calculate the 
variation of several experimentally accessible quantities 
along experimentally realizable thermodynamic paths. 
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